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ABSTRACT 

Motivated by recent observations with Herschel/PACS, and the availabihty 
of new rate coefficients for the colhsional excitation of CO (Yang et al. 2010), 
the excitation of warm astrophysical CO is revisited with the use of numerical 
i and analytic methods. For the case of an isothermal medium, results have been 

^ obtained for a wide range of gas temperatures (100 to 5000 K) and H2 densities 

(10^ — 10^ cm~^), and presented in the form of rotational diagrams, in which the 
logarithm of the column density per magnetic substate, log (Nj/gj), is plotted 
for each state, as a function of its energy, Ej. For rotational transitions in the 
wavelength range accessible to Herschel/PACS, such diagrams are nearly linear 
when n(R2) > 10® cm~^. When n(H2) ~ 10^'® — 10® cm~^, they exhibit sig- 
nificant negative curvature, whereas when n{B.2) < 10^'® cm~'^, the curvature is 
uniformly positive throughout the PACS-accessible range. Thus, the observa- 
tion of a positively-curved CO rotational diagram does not necessarily require 
the presence of multiple temperature components. Indeed, for some sources ob- 
served with Herschel/PACS, the CO rotational diagrams show a modest positive 
curvature that can be explained by a single isothermal component. Typically, 
the required physical parameters are densities in the lO'^ — 10^ cm~^ range and 
temperatures close to the maximum at which CO can survive. Other sources 
exhibit rotational diagrams with more curvature than can be accounted for by 
a single temperature component. For the case of a medium with a power- law 
distribution of gas temperatures, dN/dT oc T~^, results have been obtained for 
H2 densities 10^ - 10^ cm and power-law indices, b, in the range 1 to 5; such a 
medium can account for a CO rotational diagram that is more positively curved 
than any resulting from an isothermal medium. 



Subject headings: Molecular processes - ISM: Molecules - Infrared: ISM 
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1. Introduction 

Over the past thirty years, far-infrared line emissions from warm carbon monoxide have 
been observed from a succession of airborne and satelhte observatories: the Kuiper Airborne 
Observatory (e.g. Watson et ah 1985) the Infrared Space Observatory [ISO; e.g. Lerate et 
al. 2006), the Herschel Space Observatory (e.g. Herczeg et al. 2011) and the Stratospheric 
Observatory for Infrared Astronomy (SOFIA). Such observations have provided an invalu- 
able probe of warm molecular gas that has been heated by interstellar shock waves and 
in photodissociation regions (PDRs). In the past three years, observations performed with 
Herschel - and in particular its PACS spectrometer - have yielded a large data set of far- 
infrared CO observations, providing detections of pure rotational transitions ranging from 
J = 13 - 12 (at 200.272 /xm) to J = 49 - 48 (at 53.897 /xm). Measurements of the ladder of 
CO rotational transitions are conveniently represented by a rotational diagram, in which the 
logarithm of the column density per magnetic substate. In (Nj/gj), is plotted for each state, 
as a function of its energy, Ej. This representation is motivated by the form of the Boltz- 
mann factor, which leads - for the case of an isothermal medium in local thermodynamic 
equilibrium (LTE) at temperature, T - to a straight line for which the slope is —{1/kT); 
in this case, the y-intercept yields the total CO density divided by the partition function. 
More generally, it is convenient to define a rotational temperature, 

T,,, = -ikd\n[Nj/gj]/dEj)-\ (1) 

which - in the non-LTE case - may show a variation with Ej. 

Two features of the rotational diagrams obtained from Herschel/PACS observations 
are particularly noteworthy: (1) the plots typically exhibit positive curvature, with Trot 
increasing monotonically with Ej; and (2) the value of Trot determined from the lowest 
transitions accessible to Herschel/PACS is typically close to 300 K (e.g. van Kempen et al. 
2010a,b; Fich et al. 2010; Herczeg et al. 2011). 

The far-infrared emission spectra expected from warm CO have been modeled exten- 
sively in the literature (e.g. McKee et al. 1982, Viscuso & Chernoff 1988). The present study 
is motivated by recent observations with Herschel/PACS, as well as by the availability of 
new rate coefficients computed by Yang et al. (2011) for the coUisional excitation of CO in 
collisions with H2. I will revisit the excitation of warm CO with the use of numerical and 
analytic solutions to the equations of statistical equilibrium for the level populations. The 
calculational method is discussed in §2, including the techniques used to extrapolate the col- 
lisional rate coefficients to transitions lying higher than any considered by Yang et al. (2011). 
In §3 and §4, the predicted rotational diagrams are presented for two cases: an isothermal 
medium, and a medium in which a distribution of gas temperatures is present. Certain limit- 
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ing behaviors are accounted for with an analytic treatment, which yields excellent agreement 
with the numerical results. A brief summary follows in §5. 

2. Calculational method 

2.1. Equations of statistical equilibrium 

Using methods identical to those described in Neufeld & Melnick (1993), the equations 
of statistical equilibrium were solved to obtain the steady-state populations for 81 rotational 
states of CO (i.e. with rotational quantum number J = to 80.) As shown by Neufeld 
& Yuan (2008, Appendix), and contrary to earlier results obtained by Flower & Gusdorf 
(2009), the assumption of steady-state is generally a good one in astrophysical environments, 
because the timescale to reach equilibrium is short compared to that on which the physical 
parameters vary. Optical depth effects were treated using an escape probability method. 
The only external source of radiation was assumed to be the cosmic microwave background 
(CMB) at z = (Fixsen et al. 2009), and molecular hydrogen was assumed to be the 
dominant collision partner, with the H2 ortho-to-para ratio assumed to have its LTE value. 
With these assumptions, the level populations depend on three physical parameters: the gas 
temperature, T, the H2 density, n(B.2), and a column density parameter, A^(CO), defined by 
Neufeld & Melnick (1993). For a medium with a large velocity gradient, dv^/dz, in a single 
direction, z, N{CO) is given by N{CO)\dVz/dz\~^ . It is for this geometry that the calculation 
is actually performed, but Neufeld & Melnick (1993) have discussed how the expression for 
A^(CO) may be scaled to obtain results for different cases (including approximate results for 
media in which the large velocity gradient limit does not apply.) 

2.2. Molecular data 

I obtained the CO rotational energies from the spectroscopic constants given by Nolt 
et al. (1987), and adopted the spontaneous radiative rates obtained by Goorvitch (1994). 
Except at high densities sufficient to drive the the level populations to LTE, rate coefficients 
for the excitation of CO in collisions with H2 are of critical importance in the interpretation 
of the CO line intensities. Very recently, new estimates of these rate coefficients have been 
presented by Yang et al. (2010; hereafter YIO), who used the quantum-mechanical coupled- 
channels and coupled-states methods to compute state-to-state cross-sections, given a recent 
and highly-accurate potential energy surface obtained by Jankowski & Szalewicz (2005). 
The YIO rate coefficients covered all transitions with upper states of rotational quantum 
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number J < 40, for temperatures spanning the range 2 to 3000 K, and with both ortho- and 
para-H2 as the assumed colhsion partner. Even though this set of colhsional rate coefficients 
is the most extensive of any yet calculated, it is not sufficient to cover the large range of CO 
rotational transitions available to Herschel/PACS, and thus must be extrapolated to include 
states with J > 40. The need for extrapolation is obvious for PACS transitions higher than 
J = 40 — 39, but is also important for lower transitions, since even states with J < 40 may 
be significantly populated via a radiative cascade from states of higher J. 

I have performed an extrapolation using the general methodology adopted by Schoier et 
al. 2005; hereafter SOS) and other previous authors, in which (1) the subset of rate coefficients 
for deexcitation to the ground rotational state, 7(J(7,0), is extrapolated to states of higher 
Ju] and then (2) this subset is "expanded" using some expression based upon the infinite 
order sudden (lOS) approximation to obtain a complete set of rate coefficients, 7(Jf/, Jl)- 
However, some differences in the details of my extrapolation method are as follows: 

First, in extrapolating the set of 'j{Ju, 0), I have adopted an expression of the form 



where Eu is the energy of the upper state, Ju, and A{T) and C{T) are computed at each 
temperature from a simultaneous fit to the YIO values for Jjj = 35 and 40. This expression 
differs from that used by S05, who assumed a quadratic (rather than linear) dependence of 
ln7(J[/,0) upon J{/. A detailed examination of the results of the S05 extrapolation method 
indicate that it leads to unphysical behavior at large Ju, and thus I favored the more robust 
extrapolation expression given in equation (2). The YIO rate coefficient are shown in Figure 
1 (squares) for two example temperatures, with the extrapolation obtained from equation (2) 
plotted as dashed curves. Red and blue symbols and curves apply respectively to deexcitation 
in collisions with ortho and para-H2. The red and blue curves in Figure 2 shows the values 
obtained for C, as a function of temperature, both for ortho- (red) and para-H2 (blue). 

Second, I have used a different expression to obtain the complete set of rate coefficients, 
7( Jf/, Jl) from those obtained for J^ = 0. In the limit of large temperature, kT ^ Eu — Ei, 
Goldflam, Green & Kouri (1977, hereafter GGK) and Varshalovich & Khersonskii (1977) 
independently used the infinite order sudden (lOS) approximation to obtain the following 
relation 



7(Jc/,0) = A{T)eM-C{T){Eu/kT) 




(2) 




(3) 



where 



Ju Ju k 







) 
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is a Wigner 3-j symbol. To correct deficiencies in the lOS approximation that become 
increasingly important as kT/{Eu — Ei) decreases, McKee et al. (1982) and SOS included a 
correction factor inside the summation on the right-hand-side of equation (3), derived from 
the work of de Pristo et al. (1979). In the present work, by contrast, I have multiplied the 
entire right-hand side by a phenomenological correction factor of the form 

G{Ju,Jl) = 1+Ci[jj,j^\ [jjTj^l ■ (4) 

By optimizing the fit to rate coefficients that are available from YIO, I determined that 
the optimum values for the constants, Ci, C2 and C3 are 0.82, 2.20, and 1.22 respectively. 
Figure 3 shows the relative performance of these various methods used to "expand" the 
set of rate coefficients for cases {Jl < Ju ^ 40) in which the YIO results are available as 
a test. In each panel, the vertical axis represents the initial rotational quantum number, 
Jui and the horizontal axis shows the final quantum number, J^, with the color of each 
tile representing the ratio of the rate coefficient computed using one of the three methods 
described above to the actual value obtained by YIO: these are referred to as "GGK" (the 
use of equation 3 without correction for deficiencies in the lOS approximation), "SOS" (the 
correction introduced by de Pristo 1979 and used by McKee et al. 1982 and SOS), and 
"present work" (the correction obtained by multiplying the right-hand-side of equation 3 
by the factor G{Ju,Jl) given in equation 4). Clearly, at high temperature (T = 1000 K), 
equation (3) performs well for all transitions, but at a lower temperature (T = 100 K), the 
lOS approximation systematically underestimates the rate coefficient except when is close 
to zero or Ju. The correction used by SOS improves the fit, but the correction factor given 
in equation (4) clearly provides the best performance (even though it lacks a clear physical 
basis.) 

Figure 4 shows example rate coefficients adopted in the present work, with the vertical 
axis in each panel representing the initial rotational quantum number, Ju, the horizontal 
axis the final quantum number, Jl, and the color of each tile representing the rate coefficient. 
Tiles below the horizontal white line show the values given by YIO, while those above the 
horizontal white line show results of the extrapolation method adopted her^. 



^Thc complete set of extrapolated rate eoefScients is available from the author upon request, together 
with solutions to the equations of statistical equilibrium for the level populations; the latter were obtained 
on a fine grid of temperatures and densities, and are available as a IDL save set or a FITS file 
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3. The CO rotational diagram for an optically-thin and isothermal medium 

The optically-thin and isothermal medium represents the simplest case for which a 
rotational diagram can be determined: the fractional level populations, fj = Nj/N{CO), 
are simply computed from the equations for statistical equilibrium and then presented on a 
log- linear plot. Figure 5 shows the rotational diagrams predicted for such a medium, assumed 
in this case to be at a temperature of 1000 K. The thick curves apply to H2 densities of 10^ 
(red), 10^ (green), 10^ (blue), 10^ (magenta), and 10^ cm~^ (black), with the thinner black 
curves applying to intermediate densities spaced by 0.2 dex. The same results are presented 
in the upper panel of Figure 6, but now with the values of (fj/gj) normalized relative to 
that obtained for the lowest rotational state accessible with Herschel/PACS {viz. J = 13). 
Horizontal bars indicate the range of transitions accessible in four separate bands of the 
PACS spectrometer - denoted Long Rl, Short Rl, B2B and B2A - for the case in which the 
observations are acquired in the commonly-used "SED scan" mode (Poglitsch et al. 2010). 
The lower panel in Figure 6 shows the rotational temperature as defined in equation (1): 
Trot = -(A; din [Nj/gj]/dEj)-\ 

In the high density limit, the level populations approach LTE, the rotational diagram is 
well-approximated by a straight line, and T^ot approaches the kinetic temperature, T. As the 
density decreases, the fractional level populations for transitions in the PACS wavelength 
range drop, with the largest reductions occurring for the highest-lying transitions (which 
have the largest spontaneous radiative decay rates.) Figure 6 shows that for densities in 
the range 10®'^ to lO^cm"'^, subthermal excitation leads to a negatively-curved rotational 
diagram in which T^ot is a monotonically decreasing function of Ej. However, at yet lower 
densities, departures from LTE can introduce a positive curvature in the rotational diagram 
(i.e. the behavior that is typically observed in nature). For densities < lO^'® cm , Tj-Qi is a 
monotonically increasing function of Ej for the set of transitions accessible to PACS. 

3.1. Results in the low density limit 

Rotational diagrams for the low density limit are shown in Figure 7, with the density 
fixed at 10'^ cm~^ and for temperatures of 10^ (solid red curve), 10^'^ (green), 10'^ (blue), and 
j^gs.s (magenta). The behavior can be understood analytically by means of the following 
analysis. At low densities, most CO molecules are in the lowest two rotational states (the 
J = 1 state being populated by the absorption of CMB radiation). Any state with J > 13 
(i.e. in the range observable by PACS) is depopulated primarily as a result of spontaneous 
radiative decay, at a rate that is roughly proportional to J^, and is populated primarily 
by coUisional excitation to states with J' > J (followed by a radiative cascade except in 
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the case J' = J.) The excitation rate from J = or 1 to J' is roughly proportional to 
exp(— Xj, — Cxj'), where xj' = {Eji /kTY^'^. This last proportionality follows from the 
roughly linear behavior of the curves plotted in Figure 1. In the limit J' ^ 1, we also 
have Eji oc J'^ and gj/ oc J'. With these approximations, the equilbrium level population 
therefore varies with J in accord with 



fj 



oc 



9J Xj J^j 



xe 



-x^-Cx 



dx 



44 



e -r 



Xj—Cxj 



(2 - n'/^Ce^'^'+^/^^'eTk {xj + iC))! , (5) 



where erfc(2;) is the complementary error function and a sum over J' has been approximated 
by an integral. The corresponding rotational temperature is given by 



kTr, 



:dln [fj/gj]/dEj)-^ = 2xjEj 



Axj 



Vxj 2-7rV2Ce(^-^+^/2)'erfc(xj + iC)J 



-1 



■ (6) 



I have treated the quantity C appearing in equations (5) and (6) as an adjustable pa- 
rameter, which I have optimized for each temperature to achieve the best fit to the numerical 
results. The best-fit values for C are plotted in Figure 2 (black curve), as a function of tem- 
perature, and roug hly follo^J!| the analogous values used to extrapolate the rate coefficients 
to Jjj > 40. The dashed curves in Figure 7 present the analytic results given in equa- 
tions (5) and (6), obtained using these best-fit values of C; they evidently yield excellent 
approximations to the numerical results. 

The quantity C shows only weak dependence on temperature (Figure 2), and adopting 
a constant value of 5 yields a reasonably good approximation to all results obtained for 
T > 100 K. Because the argument of the complementary error function in equation (6) 
is always > 2.5, a further simplifying approximation is possible: for large z, eTfc{z) may 
be approximated by exp{—z'^)/{ir^^^z), the error in this approximation being < 7% for all 
z > 2.5. This additional approximation permits equation (6) to be rewritten in the form 



rot 



x^jkT 



E, 



SCxj + 2 Ej/kT + 2.5 {Ej/kry/^ + 2 ' 



(7) 



which is plotted as the dotted lines in Figure 7 (lower panel). 



^Thc black curves are systematically lower, because the curves plotted in Figure 1 show a slightly negative 
curvature: i.e. the best fit value of C for J < 20 is somewhat smaller than that ( J = 35 — 40) adopted for 
the purpose of extrapolation. 
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3.2. Comparison with PACS data 

To facilitate comparison with PACS observations, I have computed the average slopes 
in the rotational diagrams expected for each of the 4 PACS spectrometer bands: Long Rl 
(covering 140 - 220 /im), Short Rl (102 - 146 /xm), B2B (70-105 /xm), and B2A (51-73 fim). 
These average slopes lead to average rotational temperatures Tlri, Tlr2, Tb2b and Tb2A that 
may be conveniently compared with fits to the observational data. Clearly, if measurable, 
the combination of the longest and shortest wavelength spectral bands yields the maximum 
"leverage" on the physical parameters. Figure 8 shows contour plots of Tlri (red), Tb2b 
(green) and Tb2A (blue) in the logT — logn(H2) plane. While the intent of this paper is 
not to provide detailed modeling of individual sources, I note that the rotational diagrams 
observed with PACS toward embedded protostars have typically been fit as a sum of two 
components in LTE (e.g. van Kempen et al. 2010a,b; Fich et al. 2010; Herczeg et al. 2011): a 
"warm" component at Twarm ~ 250 — 400 K, and a "hot" component at Thot ~ 800 — 1400 K. 
In sources - such as NGC1333 IRAS4B (Herzceg et al. 2011) - where the observed CO 
emission extends to band B2A, and where Thot lies close to the lower end of the latter range 
(i.e. 800 K), the results in Figure 8 suggest that a single temperature component can explain 
the rotational diagram; a gas density ~ 10^ — 10^ cm~^ and a temperature ~ 3000 — 5000 K 
(i.e close to the maximum value at which CO can survive) are required to account for the 
observed rotational diagram. On the other hand, sources with larger fitted values of T^ot (e.g. 
DK Cha, observed by van Kempen et al. 2010a) have rotational diagrams with larger (more 
positive) curvature than can be accounted for by a single temperature, non-LTE component. 



3.3. Optical depth effects 

All the calculations reported in §3.1 and 3.2 above apply to optically-thin CO emission 
from an isothermal medium. CO has a dipole moment (0.11 Debye) that is much smaller 
than that of most other heteronuclear molecules, with the result that its high- J rotational 
transitions are usually optically-thin despite the large abundance of CO. Optical depth effects 
are examinecifl in Figure 9, where contour plots for Tlri (red) and Tb2A (blue) are shown for 
four values of the optical depth parameter, A(CO) (see §2.1 for the definition): 10^^, 10^®, 



■^Onc subtlety of the optically-thick case is that CO rotational diagrams are typically computed under the 
assumption that every emission line is optically-thin. In order to make a comparison with results presented 
in the literature, it is therefore necessary to obtain predictions for the level populations that would have been 
inferred under that assumption; in effect, the quantity plotted on the vertical axis is not exactly Nj/gj but, 
in reality, AirluL/lAuLgijlEij — El]), where Iul is the integrated line flux. The rotational temperatures 
plotted in Figure 9 were computed accordingly. 



- 9 - 



10^^ and 10^^ cm~^ per km/s. For a CO/H2 ratio of 10~^, non-dissociative molecular shocks 
are expected to yield ~ 10^^ cm~^ per km/s (e.g. Neufeld & Kaufman 1993); thus the 
results shown in Figure 9 suggest that optical depth effects are likely to be negligible for the 
parameters typically encountered in warm shock-heated gas. 

4. The CO rotational diagram for a continuous distribution of gas 

temperatures 

As shown above, even an isothermal medium can give rise to a CO rotational diagram 
with positive curvature - provided the gas density is sufficiently small - and may be successful 
in accounting for the CO rotational diagrams observed in some sources. Nevertheless, it is of 
interest to consider also the case where a distribution of gas temperatures is present within 
the telescope beam or along the line-of-sight. Such an admixture of temperatures might well 
occur if multiple shocks of varying velocities are present within the beam, or if shocks and 
photodissociation regions (PDRs) are both present (e.g. Visser et al. 2011). 

Evidence for multiple warm gas components of different temperatures has been inferred 
from H2 rotational diagrams observed with Spitzer, which show a positive curvature that 
cannot be accounted for by an isothermal medium (e.g. Neufeld & Yuan 2008). In addi- 
tion, the observed H2 ortho-to-para (OPR) ratio is typically considerably smaller than the 
equilibrium value of 3 expected in warm molecular gas, with a systematic tendency for the 
OPR inferred from observations of the lowest pure rotational transitions to be furthest from 
equilibrium; this behavior finds a natural explanation (e.g. Neufeld et al. 2006) if an admix- 
ture of gas temperatures is present, because the rate at which the OPR is expected to reach 
its equilibrium value is a strongly increasing function of temperature. Previous studies (e.g. 
Neufeld & Yuan 2008) have been successful in matching the observed H2 rotational diagrams 
by positing a power-law temperature distribution, dN/ dT oc T~'^db, over some range of tem- 
peratures, [Tjnin, Tlnax]- The powcr-law index, 6, typically has a best-fit value varying from 
one source to another within the range ~ 2 to 5 (Nisini et al. 2010, Yuan & Neufeld 2011, 
Giannini et al. 2011). 

In Figure 10, normalized rotational diagrams (analogous to those in Figure 6) are pre- 
sented for a power-law distribution of gas temperatures, along with the corresponding rota- 
tional temperature. These results were obtained for Tmin = 10 K (although the exact value 
of Tmin is immaterial provided it is <^ Ej/k), and for Tmax = 5000 K; my choice of Tmax is 
roughly the maximum temperature at which CO molecules can survive. The example results 
shown in Figure 10 are for a power- law index 6 = 3 and the same range of densities plotted in 
Figure 6 (and with the results for H2 densities of 10^, 10^, 10^, 10^, and 10*^ cm~^ respectively 
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appearing as thick curves in red, green, blue, magenta and black, as before). Positive curva- 
ture is predicted, both at low and high densities, with the rotational temperatures increasing 
monotonically with Ej in every case. 



4.1. Results in the high density limit 

In Figure 11, results are presented for the high-density limit. Here, the populations 
are assumed to be in LTE, and results are shown for power-law indices, b, of 1 (black), 2 
(magenta), 3 (blue), 4 (green) and 5 (red). An analytic treatment, related to that intro- 
duced by Mauersberger et al. (1988; their Appendix A) to treat the excitation of interstellar 
ammonia, can account exactly for the plotted results. For a power-law distribution of gas 
temperatures, and where the level populations are in LTE, the total population in state J 
obeys the proportionality 



9 J Jt^,„ Z^otiT) E\ 



where Z^ox.iT) oc T is the rotational partition function, zi = Ej/kT^i^^, Z2 = Ej/kT^a_^ and 
r(6, z) = J^°° t^~^e~^dt is the upper incomplete gamma function. Some further calculus and 
algebra yields the rotational temperature in the form 

Dashed black curves show the results obtained from the analytic treatment, which are in- 
distiguishable from the numerical results. For all the calculations performed here, Ej ^ 
^min = 10 K, and equation (9) becomes kT^-ot = EjT{zi,b)/T{zi,b + 1). If, in addition, 
Ej <^ Tmax = 5000 K, this expression reduces further to kT^-^t = Ej/b. The behavior in this 
limit explains the linear dependence of Trot on Ej in the left part of the lower panel. 



4.2. Comparison with PACS data 

Figure 12 is analogous to Figure 8, now showing contours of Tlri (red), Tb2b (green) 
and Tb2A (blue) as a function of the power-law index, b, and the density. By comparison with 
the isothermal case, a medium with a distribution of gas temperatures can yield rotational 
diagrams with larger positive curvature. This is manifested by larger ratios of Tb2A to Tlri, 
which provide a better fit to observational results obtained in some sources (e.g. DK Cha, 
observed by van Kempen et al. 2010a). In Figure 13, values of Tb2a/^lri are shown for the 
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two cases; the upper panel applies to a medium at a single temperature, while the lower 
panel applies to a medium with a power-law distribution of temperatures. At low density, 
an isothermal medium can result in Tb2a/^lri values as large as ~ 4, but values larger than 
4 require an admixture of gas temperatures. 

5. Summary 

1. With the use of recent collisional rate coefficients computed by YIO, and an extrapo- 
lation method discussed in §2.1, I have computed the far-infrared carbon monoxide emissions 
expected from warm astrophysical media. 

2. For the case of an isothermal medium, results have been obtained for a wide range 
of gas temperatures (100 - 5000 K) and H2 densities (10'^ — lO^cm"^), and presented 
in the form of rotational diagrams. For rotational transitions in the wavelength range 
accessible to Herschel/P ACS, such diagrams are nearly linear when n{B.2) > 10® cm~^. 
When n{B.2) ~ 10^'® — 10® cm~'^, they exhibit significant negative curvature, whereas when 
^^(1^2) < 10^'® cm~^, the curvature is uniformly positive throughout the PACS-accessible 
range. (In the intermediate density range, n{B.2) ~ 10^'® — 10^'® cm~^, there is an infiec- 
tion point where the curvature changes sign.) Thus, the observation of a positively-curved 
CO rotational diagram does not necessarily require the presence of multiple temperature 
components. 

3. For the isothermal case, an analytic treatment accounts well for the behavior in the 
low-density limit, providing simple expressions for the CO level populations and rotational 
temperatures. 

4. The average rotational temperatures expected for four PACS spectroscopic bands are 
presented as contours in the plane of density and temperature. For some sources observed 
with Herschel/P ACS, the CO rotational diagrams show a modest positive curvature that can 
be explained by a single isothermal component. Typically, the required physical parameters 
are densities in the 10^ — 10^ cm~^ range and temperatures close to the maximum at which 
CO can survive. Other sources exhibit rotational diagrams with more curvature than can 
be accounted for by a single temperature component. 

5. Optical depth effects are found to be very modest, given the typical CO column 
densities in warm astrophysical media. 

6. For the case of a medium with a power-law distribution of gas temperatures, dN/ dT oc 
, results have been obtained for H2 densities 10'^ — lO^cm"'^ and power-law indices, h, in 
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the range 1 to 5. An analytic expression is given for the CO level populations and rotational 
temperatures in the limit of high density. The average rotational temperatures expected 
for four PACS spectroscopic bands are presented as contours in the plane of density and 
temperature. 

7. A medium with a distribution of gas temperatures can account for a CO rotational 
diagram that is more positively curved than any resulting from an isothermal medium. 

It is a pleasure to acknowledge useful discussions with M. Puravankara, F. van der Tak, 
R. Visser, and D. Watson. Support for this work was provided by NASA through an award 
issued by JPL/Caltech. 

REFERENCES 

Depristo, A. E., Augustin, S. D., Ramaswamy, R., & Rabitz, H. 1979, J. Chem. Phys., 71, 
850 

Fich, M., Johnstone, D., van Kempen, T. A., et al. 2010, A&A, 518, L86 
Flower, D. R., & Gusdorf, A. 2009, MNRAS, 395, 234 
Fixsen, D. J. 2009, ApJ, 707, 916 

Giannini, T., Nisini, B., Neufeld, D., et al. 2011, ApJ, 738, 80 
Goldflam, R., Kouri, D. J., & Green, S. 1977, J. Chem. Phys., 67, 5661 
Goorvitch, D. 1994, ApJS, 95, 535 

Herczeg, G. J., Karska, A., Bruderer, S., et al. 2011. larXivillll. 07741 

Jankowski, P., & Szalewicz, K. 2005, J. Chem. Phys., 123, 104301 

Lerate, M. R., Barlow, M. J., Swinyard, B. M., et al. 2006, MNRAS, 370, 597 

Mauersberger, R., Wilson, T. L., & Henkel, C. 1988, A&A, 201, 123 

McKee, C. F., Storey, J. W. V., Watson, D. M., & Green, S. 1982, ApJ, 259, 647 

Neufeld, D. A., & Kaufman, M. J. 1993, ApJ, 418, 263 

Neufeld, D. A., & Melnick, G. J. 1991, ApJ, 368, 215 



- 13 - 

Neufeld, D. A., Melnick, G. J., Sonnentrucker, P., et al. 2006, ApJ, 649, 816 
Neufeld, D. A., & Yuan, Y. 2008, ApJ, 678, 974 
Nisini, B., Giannini, T., Neufeld, D. A., et al. 2010, ApJ, 724, 69 
Poglitsch, A., Waelkens, C, Geis, N., et al. 2010, A&A, 518, L2 

Schoier, F. L., van der Tak, F. F. S., van Dishoeck, E. F., & Black, J. H. 2005, A&A, 432, 
369 

van Kempen, T. A., Green, J. D., Evans, N. J., et al. 2010a, A&A, 518, L128 

van Kempen, T. A., Kristensen, L. E., Herczeg, G. J., et al. 2010b, A&A, 518, L121 

Varshalovich, D. A., & Khersonskii, V. K. 1977, Astrophys. Lett., 18, 167 

Viscuso, P. J., & Chernoff, D. F. 1988, ApJ, 327, 364 

Visser, R., Kristensen, L. E., Bruderer, S., et al. 2011. [aFXiv:1110.4667l 

Watson, D. M., Genzel, R., Townes, C. H., & Storey, J. W. V. 1985, ApJ, 298, 316 

Yang, B., Standi, P. C., Balakrishnan, N., & Forrey, R. C. 2010, ApJ, 718, 1062 

Yuan, Y., & Neufeld, D. A. 2011, ApJ, 726, 76 



This preprint was prepared with the AAS IATJt;X macros v5.2. 



- 14 - 



-5 



o 

cn 
O 







T = 300 K 




20 



40 

Jn 



60 



80 



o 



-10 



-15 - 



-20 - 



-25 







T = 1000 K 




20 



40 

Jn 



60 



80 



Fig. 1 - Rate coefficients for collisional deexcitation from J = Jjj ^ at temperatures of 
300 K (top panel) and 1000 K (bottom panel). Squares: results from YIO for collisions with 
ortho-H2 (red) and para-H2 (blue). Dashed lines: extrapolation adopted in the present work 
(equation 2). 
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Fig. 2 - Values adopted for C{T), when fitting the 7(J;7,0) by the expression 
A{T) exp[—C{T){Eu / kTy^'^]. Red squares: values adopted in extrapolating the rate co- 
efficients for deexcitation by ortho-H2. Blue squares: values adopted in extrapolating the 
rate coefficients for deexcitation by para-H2. Black squares: values that yield the best fit to 
the rotational diagram (as approximated by equation (5) in §3.1). 
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Fig. 3 - Relative performance of various methods (see text) used to "expand" the set of rate 
coefficients for cases (J^ < Ju ^ 40) in which the YIO results are available as a test. In 
each panel, the vertical axis represents the initial rotational quantum number, Ju, and the 
horizontal axis shows the final quantum number, Jl. The color of each tile indicates the 
ratio of the rate coefficient computed using each method to the actual value obtained by 
YIO. 
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Fig. 4 - Rate coefficients for collisional deexcitation of CO by H2, with the vertical axis in 
each panel representing the initial rotational quantum number, J^/, the horizontal axis the 
final quantum number, J^, and the color of each tile representing the rate coefficient. Tiles 
below the horizontal white line: values given by YIO. Tiles above the horizontal white line: 
extrapolation described in §2.1 
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Rotational diagram at 1000 K 
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Fig. 5 - Rotational diagrams predicted for an optically-thin, isothermal medium at a temper- 
ature of 1000 K. The thick curves apply to H2 densities of 10^ (red), 10^ (green), 10^ (blue), 
10^ (magenta), and 10^ cm~^ (black), with the thinner black curves applying to intermediate 
densities spaced by 0.2 dex. 
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Fig. 6 - Upper panel: same as Figure 5, but with the level populations normalized 
to those for J = 13. Bottom panel: corresponding rotational temperatures, Trot = 
— (/cdln [Nj/ gj]/dEj)~^ . Horizontal red and blue bars indicate the range of rotational levels 
probed by observations in each of the four labeled PACS spectral bands. 
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Normalized rotational diagram at n(H2) - lO^cm' 
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Fig. 7 - Rotational diagrams (top panel) and rotational temperatures (bottom panel) pre- 
dicted for an optically-thin, isothermal medium at an H2 density of lO'^cm"'^, with solid 
red, green, blue and magenta curves applying to temperatures of 10^, 10^'^, 10'^ and 10^'^ K 
respectively, and the approximate analytic treatment (equation 5) shown as dashed curves. 
The dotted curves were obtained with a further simplifying approximation (equation 7; see 
text). 
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Fig. 8 - Contour plots of the average rotational temperatures (in Kelvin) predicted for CO 
transitions in the Long Rl (140 - 220 /im), B2B (70-105 /im), and B2A (51-73 /im) bands 
of Herschel/P ACS, as a function of temperature and H2 density. These are shown in red 
(Tlri), green (Tb2b), and blue Tb2A, respectively. Results apply in the optically-thin limit. 
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Fig. 9 - Same as Figure 8, but for four different values of the optical depth parameter, 
iV(CO) (see text). 
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Normalized rotational diagram for b = 3.0 
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Fig. 10 - Rotational diagrams, and rotational temperatures, Trot = — (fcdln [Nj/gj]/dEj)~^, 
predicted for an optically-thin medium with a power-law distribution of temperatures of the 
form dN/dT = T~^dT over the range 10 to 5000 K. Thick curves apply to H2 densities of 10^ 
(red), 10^ (green), 10^ (blue), 10^ (magenta), and 10^ cm~^ (black), with the thinner black 
curves applying to intermediate densities spaced by 0.2 dex. The power-law index assumed 
here is h = 3.0, and the level populations in the top panel are normalized with respect to 
those for J = 13 
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Normalized rotational diagram for n(H2) — 10^ cm' 
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Fig. 11 - Rotational diagrams (top panel) and rotational temperatures (bottom panel) pre- 
dicted for an optically-thin medium, with a power-law distribution of gas temperatures, at an 
H2 density of 10^ cm~^. Solid black, red, green, blue and magenta curves apply to power-law 
indices, b (see text), of 1, 2, 3, 4, and 5 respectively. Results from an approximate analytic 
treatment (equation 9) are shown as dashed curves. 
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Fig. 12 - Contour plots of the average rotational temperatures (in Kelvin) predicted for 
CO transitions in the Long Rl (140 - 220 /xm), B2B (70-105 /xm), and B2A (51-73 /im) 
bands of Herschel/PACS, as a function of the power-law index, b, and the H2 density. These 
are shown in red (Tlri), green (Tb2b)5 and blue (Tb2a)5 respectively. Results apply to an 
optically-thin medium with a power-law distribution of gas temperatures. 
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Fig. 13 - Red contours: average rotational temperatures (in Kelvin) predicted for CO transi- 
tions in the PACS Long Rl band (red). Black contours: ratio of average rotational tempera- 
ture in the PACS B2A band to that in the Long Rl band. Results apply to an optically-thin 
medium at a single temperature (top panel) or with a power-law distribution (see text) of 
gas temperatures (bottom panel). 



